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ABSTRACT 

Spicules were recently found to exist as two different types when a new class of so-called type 
II spicules was discovered at the solar limb with the Solar Optical Telescope onboard the Hinodc 
spacecraft. These type II spicules have been linked with on-disk observations of Rapid Blue-shifted 
Excursions (RBEs) in the Ha and Ca II 8542 lines. Here we analyze observations optimized for the 
detection of RBEs in both Ha and Ca II 8542 lines simultaneously at a high temporal cadence taken 
with the CRisp Imaging SpectroPolarimeter (CRISP) at the Swedish Solar Telescope (SST) on La 
Palma. In this study we used a high-quality time sequence for RBEs at different blue-shifts and 
employed an automated detection routine to detect a large number of RBEs in order to expand on 
the statistics of RBEs. We find that the number of detected RBEs is strongly dependent on the 
associated Doppler velocity of the images on which the search is performed. Automatic detection of 
RBEs at lower velocities increases the estimated number of RBEs to the same order of magnitude 
expected from limb spicules. This shows that RBEs and type II spicules are indeed exponents of the 
same phenomenon. Furthermore, we provide solid evidence that Ca II 8542 RBEs are connected to 
Ha RBEs and are located closer to the network regions with the Ha RBEs being a continuation of 
the Ca II 8542 RBEs. Our results show that RBEs have an average lifetime of 83.9 seconds when 
observed in both spectral lines and that the Doppler velocities of RBEs range from 10 to 25 km s _1 in 
Ca II 8542 and 30 to 50 km s" 1 in Ha. In addition, we automatically determine the transverse motion 
of a much larger sample of RBEs than previous studies, and find that just like type II spicules, RBEs 
undergo significant transverse motions of order 5-10 km s _1 . Finally, we find that the intergranular 
jets discovered in Big Bear Solar Observatory are a subset of RBEs. 



1. INTRODUCTION 

Observing the chromosphere at the solar limb, a great 
number of thin jet-like features called spicules can be seen 
protruding outwards from the limb. While spicules have 
been known a nd studied for a long time (for an extensive 
overview, see iBeck ers 1968; Sterling 2000), progress in 
understanding their nature has been relatively modest. 
However, interest for this phenomenon has revived con- 
siderably in recent years. The seeing-free environment 
and aperture of sufficient size for high spatial resolution 
allowed the Hinode satellite to provide an unprecedented 
view of the solar limb in the form of high temporal res- 
olution, long duration time series of Ca II H filtergrams. 
From th ese observations, a new class of spicules was dis- 
covered (|De Pontieu et al. 2007) that is characterized by 
short life times (typically 10-100 s), vigorous dynamics, 
and exclusively upward motion (50-150 km s" 1 ). Their 
nature is quite different from what are now considered 
as classical spicules, or type I spicules, which have life 
times from 5 to 10 min, and display clear rise and down- 
fall phases at more modest velocities (20-30 km s _1 ). 
The characteristic sideways swaying motion that this 
newly identified class (or type II spicules) displays, has 
been regarded as a sign that the chromosphere is per- 
meated with Alfvenic waves of sufficient energy to ac- 
celerate _ttie_jc4a£_wajidand_potentially heat the quiet 
corona (|De Pontieu et al.ll2007| ). Later studies have fur- 
ther strengthened the conjecture that type II spicules 
have an important role in mass loading and heating of 



the corona (|De Pontieu et al.ll2009l MHh 

The driving mechanism behind type I spicules seems 
to be well understood. Type I spicules have been 
identified as the off-limb manifestation of what is 
known as dynamic fibrils in active re gions on the disk 
The o b servations and modeli ng of 
(|205l . lHansteen etHl (ppM ) and 
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(2007) demonstrated that dynamic fibrils are driven by 
magnetoacoustic shocks that result from flows and waves 
that leak along magnetic field lines from the photo- 
sphere into the chromosphere. A similar mechanism is 
responsible for driving part of th e quiet Sun mottles 
(jRouppe van der Voort et aT1l2007h which helps explain 
why type I spicules can be detected in both active regions 
and quiet Sun. 

In contrast to type I spicules, the formation and driver 
of type II sp icules is not well under s tood and still un- 
der debate. iMartmez-Svkora et al.l (|201 lh describe in 
detail the formation of what resembles a type II spicule 
in their 3-dimensional radiative magneto-hydrodynamics 
simulation. The simulated jets shows similarities to type 
II spicules in that rapid ejection of cool chromospheric 
plasma into the corona occurs while it is being heated 
to temperatures that are high enough to cause apparent 
fading in typical chromospheric diagnostics. In the simu- 
lation, small-scale flux emergence leads to chromospheric 
plasma being accelerated by a strong, mostly horizontal, 
Lorentz force into a region of very strong vertical mag- 
netic field. The subsequent increase in gas pressure in- 
duces the ejection of plasma into the low-density corona 
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while being heated to coronal temperatures by heating 
processes at the footpoint. 

While this scenario is tantalizing, further modeling is 
required in order to firmly establish the physics behind 
type II spicules. Perhaps more importantly, high-quality 
observations are essential in order to constrain the mod- 
els and guide the direction of investigation. A compli- 
cating factor here is that due to their narrow spatial ex- 
tents, short lifetimes and significant displacement, as well 
as the superposition of many spicules along the line of 
sight, spicules are notoriously difficult to properly mea- 
sure at the limb. Finding the disk-counterparts of type 
II spicules makes it possible to separate them spatially 
and overcome the problem of superposition encountered 
a t the limb. 

iLangangen et al.l (|2008f ) investigated so called "rapid 
blueshifted excursions" (RBEs) in on-disk Ca II 8542 
data obtained with the Interferometric BIdimensional 
Spectrometer (IBIS) at the Dunn Solar Telescope. In 
their data, events with short-lived blue-ward asymme- 
tries of the Ca II 8542 line were tentatively linke d to 
type II spicules. iRouppe van der Voort etall (|2009L Pa- 
per I) used high-quality Ha and Ca II 8542 observations 
from the Crisp Imaging Spectropolarimeter (CRISP) at 
the Swedish Solar Telescope (SST) to make a firm con- 
nection between these RBEs and type II spicules. They 
found a large number of RBEs for which properties such 
as lengths, velocities, and life times were found to agree 
well with what was measured for type II spicules. In ad- 
dition, the appearance, including both swaying motions 
and apparent propagation away from magnetic field con- 
centrations in the photosphere, was considered as a clear 
indication that RBEs are the disk counterpart of type II 
spicules. 

Recently, Uudge et al.l fj201 lh challenged the identifi- 
cation of RBEs as type II spicules by arguing that the 
global occurrence rate of RBEs as inferred from the Pa- 
per I measurements is too low. In this work, we further 
investigate the occurrence rate of RBEs in new high- 
quality CRISP observations as well as the observations 
used in Paper I. For the new observations, we employed 
a dual-line program that provides both spectrally well 
sampled Ca II 8542 and Ha profiles at relatively high 
temporal cadence. This program allows for a direct com- 
parison between RBEs measured in Ca II 8542 and Ha, 
a comparison that could only be done indirectly from 
the sequentially recorded data sets used in Paper I. Fur- 
thermore, we expand on previous results (from Paper I) 
by extending the automated detection routine for RBEs 
so that we can automatically determine life times and 
transversal motion for a much larger statistical sample 
than in Paper I. 

2. OBSERVATIONS AND DATA REDUCTION 

The data was obtained from obser vations at the 
Swedis h 1-m Solar Telescope (SST, Scharmc r et al.l 
12003a!) on La Palma us ing the CRisp I magin g Spec- 
troPolarimeter (CRISP, IScharmer etall |2008| ) instru- 
ment. CRISP is equipped with a dual Fabry-Perot inter- 
ferometer system and three high-speed low-noise CCD 
cameras. These cameras operate at a frame rate of 35 
frames per second with an exposure time of 17 ms per 
frame and are synchronised by an optical chopper. Two 
of the cameras are positioned behind the FPI and a po- 



larizing beam splitter, while the third camera, which is 
used as an anchor channel for image processing, is po- 
sitioned before the FPI but after the CRISP pre-filter 
and is referred to as the wide-band channel. The field 
of view (FOV) with the CRISP instrument is roughly 
61x61 arcseconds with an image scale of O'.'0592/pixel. 
CRISP allows for fast wavelength tuning (<50 ms) be- 
tween any two positions within a spectral region given by 
the spectral width of the prefilter, which makes it ideal 
for studies of the dynamic time evolution of the chro- 
mosphere through imaging spectroscopy. Here we are 
interested in the Ha and Ca II 8542 spectral lines. For 
Ha, CRISP has a transmission FWHM of 66 mA and 
pre-filter FWHM of 4.9 A. For Ca II 8542, the transmis- 
sion FWHM is 111 mA and pre-filter FWHM is 9.3A. 
High spatial resolution down to the telescope diffraction 
limit (\/D=tf.'U and C/18 for Ha and Ca II 8542 respec- 
tively) is accomp lished by use of the SS T adaptive op- 
tics system (SST, IScharmer et al.ll2003bl) and the Multi- 
Object, Multi-Frame Blind Deconvolution image r estora- 
tion technique (MOMFBD lvan Noort et al.ll2005l ). 

A 38-minute dataset obtained on June 27, 2010, (17:35- 
18:13 UT) has been analyzed. The target area was a 
coronal hole with a unipolar magnetic island, at solar 
coordinates (x,y)«(-70",508"), (^=0.84). Figure[T]gives 
an overview of the target region. 

The observational program was made specifically for 
the detection of RBEs in Ca II 8542 and Ha, simultane- 
ously, with a rather high cadence so that the time evo- 
lution of RBEs could be resolved properly. The prefilter 
wheel limited the amount of dead time to 0.63 s when 
switching between the two lines. For the Ca II 8542 
line, we ran a program with symmetric sampling at 7 
positions in each of the two wings, from ±120 to ±1452 
with steps of 222 mA. In addition, we observed the line 
core and the blue wing at —600 mA with modulation of 
the liquid crystals with the purpose of obtaining photo- 
spheric magnetograms. For Ha we ran a program with 
dense sampling of 28 positions in the blue wing, from 
-300 mA to -2082 mA with steps of 66 mA. More de- 
tailed information on the sampling positions and their 
respective Doppler velocities can be seen in Fig. [2] d and 
c, respectively. This program was based on our experi- 
ences from Paper I and specifically designed to simulta- 
neously detect RBEs in both lines. RBEs are detected 
in Ca 8542 in Dopplermaps constructed from the sub- 
traction of symmetric positions from line centre, and in 
Halpha in maps constructed from the subtraction of the 
far blue wing image from positions closer to line centre. 
The program is sufficiently fast (11.8 s cadence) to allow 
for temporal analysis of RBEs. The symmetric sampling 
of the Ca 8542 line allows for verification that the de- 
tected events are indeed asymmetric in the blue wing 
which is the defining characteristic of RBEs. 

We acquired 8 exposures per line position which were 
used for Multi-Object, Multi- Fra me Blind Deconvolutio n 
image restoration (MOMFBD Ivan Noort et all 120051) . 
Exposures from each spectral line were processed in sepa- 
rate MOMFBD restorations. Precise alignment between 
the wideband and narrowband cameras is achieved by a 
separate alignment procedure involving a reference pin- 
hole array target. For the MOMFBD restoration of the 
sequentially recorded CRISP images, the wideband chan- 
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Figure 1. Overview of the observed target area, (a): Halpha line core image, (b): Ca II 8542 line core image with overplotted in rainbow 
color scale a density map of identified RBEs detected at —20 km s^ 1 over the whole time series, (c): Ha wideband image (FWHM=4.9A 
CRISP prefilter) with overdrawn temporally resolved RBEs that are identified in both Ca II 8542 and Ha: the Ca II 8542 RBEs are drawn 
in red, the Ha RBEs are green, the region where they overlap is yellow, (d): Ca II 8542 —600 mA Stokes V image in gray color scale with 
overplotted a density map of the number of identified RBEs detected in Ha at —45 km s~ 1 during the whole time series. 



nel serves as a so-called anchor channel that ensures ac- 
curate alignment between the different line positions. For 
certain seeing conditions, the assumption of the size of 
the iso-planatic patch used to divide the FOV in subim- 
ages (here 128 x 128 pixels, or 7'.'6 squared) is not suf- 
ficiently accurate and the derived PSF is not the opti- 
mal solution over the whole area of the subimage. This 
results in remaining rubber sheet deformations between 
the CRISP line positions in a single restoration. To re- 



duce this effect, we employed an extra step in the pro- 
cessing following an idea from Henriques (private com- 
munication). The wideband images are used twice in a 
restoration: 1. all images combined in the normal way 
resulting in a single restored "anchor" image, 2. sepa- 
rated in sets associated with the corresponding CRISP 
line positions resulting in one restored wideband image 
per CRISP line position. Even though individual wide- 
band exposures enter twice in the restoration procedure, 
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they are only used once for the determination of the 
PSF. For the use of the wideband images separated in 
sets per CRISP line position, the weight is put to zero 
for the wavefront sensing. The separately restored wide- 
band images have the same rubber sheet deformations as 
the corresponding CRISP line positions which can now 
be accurately measured by cross-correlation to the an- 
chor wideband images. The deformation grid determined 
from the wideband images is then applied to the CRISP 
images which results in precise alignment between the 
different line positions. The data analyzed in this study 
benefited significantly from this extra step in the pro- 
cessing. 

After the MOMFBD restoration the images from the 
individual line scans were combined to form a time series 
and the Ca II 8542 images in each scan were aligned 
with respect to the Ha images from the same time step. 
This alignment between the two lines was done using 
the wideband images which show photospheric scenes for 
both spectral lines so that the alignment is done with 
high accuracy. The images from both lines were then 
de-rotated to account for diurnal field rotation, aligned 
and dc-strctched to remove warping due to seeing effects 
by determining local offsets on the wideband images and 
applying these to the CRISP images. 

For the Ca II 8542 —600 mA magnetograms, 
all 32 exposures of the four different liquid crys- 
tal states were included in a single MOMFBD 

described by 
and 



restoration , following the scheme 



van Noort fe Rouppe van der Voortl 
Schnerr et al.l (|2011| ). For the polarimetric response 
of the teles cope at 8542 A, w e used the model 
derived bv Ide la Cruz Rodri'guezI (120101) (also see 



Ide la Cruz Rodriguez fc Socas-Navarroll20lil ). With the 
17 ms exposure time, the effective acquisition time for 
the resulting Stokes V magnetogram is 544 ms. The 
noise level is estimated at 0.14% of the continuum. 
It should be noted that at -600 mA, the Ca 8542 
magnetogram is dominated by photospheric signals and 
is not affected by strong Doppler shifts. Strong Doppler 
shifts make single wavelength magnetograms closer to 
the line core, that are more chromospheric in nature, 
difficult to interpret. Furthermore, it should be noted 
that the Ca II 8542 line is less sensitive to polarization 
than for example more traditional diagnostics, like the 
Fe I 630 lines. In addition, the wavelength position used 
here, —600 mA, is relatively less sensitive to polarization 
than positions closer to the line core. All these factors 
combine to a rather limited diagnostic value of the 
magnetograms used in this study, in particular for the 
weak magnetic field environment of the coronal hole 
observed here: the strongest signals barely reach 1% po- 
larization. Despite these limitations, the magnetograms 
serve as a reliable indicator of the position and polarity 
of the strongest magnetic field concentrations. These 
concentrations appear as photospheric bright points 
in continua and in the far wings of the spectral lines 
studied here. 

We acquired another high-quality data set with the 
same observing program during the 2010 SST observing 
campaign: on 03-Jul-2010, we observed a coronal hole at 
(x,y)«(19,862) (/i=0.41) for 45 min. In this paper, we 
present results from the analysis of the 27-Jun-2010 data 



set but the 03-Jul-2010 data set is used as a reference to 
verify and check our results and conclusions. 

3. METHOD 

Ha Co II 854.2 nm 
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Figure 2. Details of an RBE observed in both Ha (left column) 
and Ca II 8542 (right column). Panels (a) and (b) show the tem- 
poral evolution of the spectrum at a central location in the RBE, 
marked with asterices in panels (e) and (f) (note that this is not the 
same spatial location for the two spectral lines). The temporal axis 
is centred on the time used for the panels below. Panels (c) and 
(d) show the spectral profile at the asterix' location with the solid 
line, the average spectrum over the full FOV with the solid line 
with diamonds at the sampling points, and the difference between 
these profiles as the dashed line. Panels (d) and (e) show blue wing 
images at -45 km s _1 for Ha and —20 km s _1 for Ca II 8542. The 
RBE is seen as a dark streak which was seen to originate from the 
magnetically enhanced region in the upper right direction in the 
time steps before. The Ca II 8542 RBE is shorter and closer to 
this region than the Ha RBE. Panels (g) and (h) show the maps 
used for automated detection: a difference image in Ha with the 
—45 km s _1 image subtracted from the far wing —2 A image, and 
a Ca II 8542 Doppler map at ±20 km s" 1 . The detection skeleton 
from the automated method is marked with a thin black line. 



The dataset was searched for RBEs in the Ca II 8542 
and Ha lines separately, where they appear in the spec- 
tral wings as roundish or elongated dark absorption fea- 
tures against the photospheric background. An auto- 
mated algorithm that locates RBEs by isolating the loca- 
tion and time of a long thin feature at high blueshifts was 
used to find a large number of RBEs in both datasets and 
both spectral lines in order to get reliable RBE statistics. 
Apart from some minor refinements, the method is iden- 
tical to the method used for Paper I. 

To make a clean detection of the upper chromospheric 



5 



RBEs, the background, which often contains linear pho- 
tospheric features (such as intergranular lanes), has to 
be reduced or removed from the images. For Ca II 8542 
this is done by subtraction of the opposite wing position, 
thereby constructing a Doppler map, and for Ha this is 
done by subtracting the far blue wing image, at 2A, cre- 
ating a difference image. To minimise the number of false 
detections, a lower limit for the length of blue-shifted 
structures is set at ^725 km, 17 pixels, and any thin 
blue-shifted feature with a length below this threshold is 
not counted. Like in Paper I, we find that this length 
threshold limits the probability of false detections signif- 
icantly. A shorter length threshold allows small linear 
features in intergranular lanes to enter the sample which 
upon further scrutiny do not display the typical spectral 
character of RBEs. The output of the automated algo- 
rithm contains every detected event in every timeframe, 
and it also attempts to connect these single events over 
several timeframes by making chains of events which rep- 
resents the time evolution of an RBE over several time 
steps. If an event from one timeframe has 30% of its pix- 
els within a 3 pixel distance from any pixel in an event 
from the next timeframe, these two events are considered 
to be the same RBE observed at different stages in its 
time evolution. Because the lengths of events may vary 
greatly from frame to frame, the check has to be done 
in both temporal directions as the following event, which 
is supposed to link up with the current event, might be 
shorter than 30% of the current event's length. In this 
case, even if the following event lies completely on top of 
the current event, it would not be recognised as a valid 
link, while when performing the check " backwards" there 
is a 100% overlap and the link is approved. This check 
is performed on all events in a timeframe using all the 
events from the next timeframe and also the frame af- 
ter that, meaning it allows for a jump in the chain of 
one timestep thus accounting for occasional moments of 
bad seeing. This automated method of connecting RBEs 
in time is an improvement compared to Paper I, where 
the lifetimes of RBEs were estimated manually by look- 
ing at a subset of 35 separate events and following their 
existence in following time steps. 

Having the chains made automatically gives a large 
sample from which one can study the time evolution and 
make accurate statistics of the typical lifetimes of RBEs. 
The cases where RBEs are connected in chains are also 
less likely to be false detections, as it would take errors to 
occur very close together spatially and temporally. When 
it comes to the selection threshold for determining the 
temporal connection in chains, the threshold values were 
chosen based on trial and error. Allowing for events to be 
further away from one another introduced a significant 
number of wrong connections, particularly in the regions 
where the RBE detections are densest, and by constrain- 
ing the percentage of pixels that has to be within the 
selected distance, too many good connections were dis- 
allowed. We found that the chosen parameters resulted 
in the most reliable statistics. 

Once the chains are found, these can be used to get 
a measure of the transverse displacement and velocity. 
This is achieved by finding the average orientation of all 
events in a chain and determine its perpendicular direc- 
tion. The transverse distance is then measured by tak- 



ing the middle point of each event and projecting it onto 
the normal which allows for both a measurement of the 
distance traversed transversally from one event to the 
next and also the velocity at which the RBE is moving 
transversally as we know the time difference between the 
two events. As explained above, there is an upper limit 
of 3 pixels on the spatial separation between RBEs in 
subsequent time frames. Combined with the cadence of 
11.8 s, one would expect to have an upper limit for the 
transverse velocity of 11 km s . However, this is not the 
case since the transversal displacement is measured from 
the RBE midpoints along the normal to the average di- 
rection of the RBE and this can give a larger separation 
than 3 pixels. 

We have made extensive use o f CRISPEX 
(jVissers fc Rouppe van der Voortl 120121 ) a widget 
based analysis tool programmed in the Interactive Data 
Language (IDL) that allows for efficient exploration 
of multi-dimensional datasets. With CRISPEX, the 
Ha and Ca II 8542 datasets can be examined at the 
same time and the tool was particularly useful for 
the verification of the RBE identifications from the 
automated detection method. 

Figure [2] shows details of one RBE event, clearly iden- 
tified in both Ha and Ca II 8542. Panels (e) and (f) show 
the RBE as dark streaks against a background of granu- 
lation in Ha and inverse granulation in Ca II 8542. For 
the spatial locations marked with asterisks, the spectral 
time evolution is shown in the top panels. This par- 
ticular RBE displayed little sideways motion so that a 
single pixel location shows the actual spectral evolution 
of the RBE. These panels illustrate the defining charac- 
teristic of RBEs: a short-lived asymmetry in the blue 
wings of these strong absorption lines. The bottom pan- 
els show the maps used in the automated detection pro- 
cess: in Ha the disturbing background is attenuated by 
subtraction of the corresponding far blue wing image. 
Through this subtraction, the granulation intensity pat- 
tern is inverted so that intergranular lanes, narrow linear 
features, become bright. Our detection method is based 
on the negative of this difference images and searches for 
bright linear structures, removing the risk of intergranu- 
lar lanes becoming false detection as these are now dark. 
In Ca II 8542, a Dopplermap is used and RBEs are iden- 
tified as bright linear structures against a low-contrast 
background. The detected structures are marked as thin 
black lines. Panels (c) and (d) show the detailed spec- 
tra at the time and location of the asterisks in panels 
(e) and (f) as solid lines. To separate the RBE spectral 
profile from the background profile, a reference spectrum 
constructed from averaging over the whole FOV is sub- 
tracted, resulting in the dashed profile. This RBE profile 
is clearly a blue-shifted component separated from the 
background spectrum, a characteristic that is also ap- 
parent from its dynamical evolution which appears to be 
unrelated to the chromospheric background seen closer 
to line centre. The RBE profile is used to determine 
the Doppler velocity and Doppler width, following the 
method of Paper I. For this location in the RBE, the 
Ha Doppler velocity is 39 km s _1 and width 16 km s -1 , 
the Ca II 8542 Doppler velocity is 23 km s _1 and width 
13 km s _1 . 

The length of an RBE was determined using the same 
method as Paper I: by first extending the detection skele- 
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Figure 3. Histograms of the number of detected events as a func- 
tion of the Doppler velocity used in the detection image. Ho detec- 
tion numbers in solid lines, Ca II 8542 in dashed lines. The right 
axis scale indicates the estimated total number of events at any 
time over the whole surface of the Sun. 

ton (the thin black line in panels (g) and (h) outlining 
the detected event) by 20 pixels in both directions and 
dilation in width and then determining the points where 
either the Doppler velocity or Doppler width of the RBE 
profile is 0. This definition is equivalent to measuring 
the length for which a feature shows enhanced absorp- 
tion compared to the average spectral profile. Hence, 
detected events have a possibility of becoming shorter 
than the minimum length threshold set by the detection 
routine since the measured Doppler width can become 
even though there is signal in the detection map. 

4. RESULTS 
4.1. Number of RBE detections 

We found that the number of RBE detections varies 
significantly with the Doppler velocity from which the de- 
tection maps are constructed. Figure [3] shows histograms 
for the number of detected events as a function of the 
Doppler velocity in the detection map. For both spectral 
lines, we see a significant increase in the number of de- 
tected events for lower Doppler velocities. This trend is 
most pronounced for Ha, where we measure an increase 
from 341 at 63 km s^ 1 to 9285 events at 30 km s — 1 . 
We manually verified a large subsample of the events de- 
tected at 30 km s _1 and we confirmed that almost all 
events are genuine RBEs. For Ca II 8542 we see a raise 
in the number of detected events from 346 at 35 km s _1 
to 4116 at 12 km s _1 . Manual inspection at 12 km s _1 
however revealed a high fraction of false identifications. 
We conclude that this line position cannot be trusted for 
the automated detection of RBEs. For the remainder of 
this paper, we present results from Ca II 8542 RBEs de- 
tected at 20 km s" 1 , where we find the highest number of 
detections (2717 events). For Ha, we present results from 
RBEs detected at 45 km s _1 (2099 events) which gives 
a sample of similar size as for Ca II 8542. This choice 
is rather arbitrary since we regard all Doppler velocities 
reliable for RBE detection. We found it practical to have 
similar sample sizes for the two lines for the process of 
linking up the Ha and Ca II 8542 RBEs. 

The right axis scaling of Fig. [3] gives an estimate of 



the total number of RBEs at any time on the Sun. This 
estimate is derived from the average number of detections 
per time step extrapolated to the whole surface of the 
Sun. 

4.2. Statistical properties of RBEs 

Figure [Hi shows a density map of the spatial locations 
of the 2099 Ha RBEs detected at 45 km s _1 drawn on 
a Stokes V magnetogram. Most RBEs are found around 
the patch of unipolar magnetic field in the upper right 
part of the FOV, the region that is dominated by fibril 
structures in the line core images and bright points in the 
wideband image. There is another notable concentration 
of RBE detections below the centre which is associated 
with a small concentration of enhanced stokes V signal, 
and another smaller concentration to the left edge of the 
FOV. Besides these RBE concentrations there are a num- 
ber of isolated RBEs scattered over the FOV that are not 
obviously associated with patches of stokes V signal. As 
explained in Sect. 2, the diagnostic potential of the Ca II 
8542 -600 mA magnetograms is rather limited and it is 
well conceivable that a more sensitive polarization diag- 
nostic would reveal association of the more isolated RBE 
signals with areas of enhanced polarization. 

Of the 2099 Ha events, 498 were detected only in one 
timeframe, while the rest linked up to other events in 
nearby timeframes to create a total of 453 RBEs existing 
over several sequential time steps. We will call the events 
that could not be linked up with events in other time 
steps "single-frame RBEs" , and the others "multi-frame 
RBEs" . 

The map of Ca II 8542 events (Fig. [TJd) shows a simi- 
lar distribution of most RBEs being associated with the 
magnetic field concentrations and a number of events 
scattered over the the FOV. Of the 2717 Ca II 8542 
events, 827 were single-frame events, and the rest was 
linked up as 507 temporally resolved or multi-frame 
RBEs. 

Most of the multi-frame and single-frame RBEs were 
located in regions with enhanced polarization signal, and 
of the isolated RBEs that lie scattered over the FOV, the 
majority are single-frame RBEs, as can be derived from 
Fig. [Tb and d. 

For each individual detected event the Doppler veloc- 
ities and Doppler widths were calculated at each pixel 
in an event. The position of the pixels and timeframe 
in which it is found were used for a measure of lengths 
and lifetimes of the RBEs. Figure |4] and [5] display his- 
tograms of all the important properties of the 1890 and 
1601 events that were found to make up 507 and 453 
multi- frame RBEs in Ca II 8542 and Ha, respectively. 
The single-frame events are not included in this statistic 
as they lack time evolution, and thus cannot contribute 
to some of the histograms, such as transverse motions. 

In panel (a) of Fig.@]thc lengths of all events are shown 
for both Ca II 8542 and Ha. Typically, the Ha RBEs 
are found to have lengths 2—6 Mm, with an average of 
3.5 Mm, while the Ca II 8542 RBE lengths are seen to lie 
between 1.5 and 4.5 Mm with an average of 2.9 Mm. A 
few Ca II 8542 RBEs were found to have lengths shorter 
than the detection threshold of 725 km which is possible 
since the method of measuring the RBE length is dif- 
ferent from the RBE detection (see Sect. [3]). No short 
Ha RBEs were found and we see a sharp lower cutoff at 




Figure 4. Histograms for (a) length, (b) Doppler velocity, (c) Doppler width, and (d) lifetime of RBEs that were identified for more than 
1 time step. Solid lines are He* RBEs, dashed lines Ca II 8542 RBEs. The histograms for Doppler velocity and width are average measures 
for each time step in each event. 



2 Mm. 

Panel (b) of Fig. 0] shows the average Doppler velocity 
of all RBEs. The average RBE Doppler velocities lie 
between 10-25 km s" 1 for Ca II 8542 and 30-50 km s" 1 
for Ha, with the mean value of the distribution of average 
Doppler velocities at 17.7 and 36.8 km s _1 , respectively. 
The maximum of all Doppler velocity measurements is 
as high as 50 km s _1 for Ca II 8542 and 62 km s _1 for 
Ha. 

Panel (c) of Fig. [4] shows the average Doppler width 
of all RBEs. As for the Doppler velocity, the Doppler 
widths are noticeably different in Ca II 8542 and Ha with 
averages of 10 km s~ 1 and 15.3 km s _1 , respectively, and 
the range going from to 17 km s _1 in Ca II 8542, while 
Ha RBEs display Doppler widths from 7 to 23 km s _1 . 

Panel (d) of Fig. [4] shows the lifetimes of the RBEs. 
There is no major difference between the lifetime of 
Ca II 8542 and Ha RBEs and they both have lifetimes 
varying from 24 s up to 2.5 min. There were a few detec- 
tions of very long lived RBEs with lifetimes of more than 

3 min. They were all found in regions with a high density 
of RBEs and visual inspection exposed these events as 
falsely connected chains of several separate multi-frame 
RBEs. 

Figure [5] show different measures of the transverse mo- 
tion of RBEs. Panel (a) shows a histogram of the max- 
imum transverse displacement which is the distance be- 
tween the two most extreme transverse positions of the 



RBE during its lifetime. All RBEs have some transverse 
displacement and several of the RBEs display transverse 
displacements of more than 400 km. The scatter plot in 
panel (b) shows the integrated transverse displacement, 
or the sum of the absolute values of each displacement 
during the RBE lifetime, against the maximum trans- 
verse displacement. The dotted line marks where the 
integrated transverse displacement is 50% larger than 
the maximum transverse displacement, the RBEs that 
lie above this line display considerable back-and-forth 
swaying motion during their lifetime. The RBEs that 
lie close to the solid line are mainly moving in one direc- 
tion. Panel (c) shows the maximum transverse velocity 
for each RBE and panel (d) shows all measured trans- 
verse velocities. Transverse motions can be at velocities 
as high as 20 km s" 1 although most events' velocities 
lie below 15 km s _1 with the bulk of the peak velocities 
around 5 km s _1 . 

We observe many RBEs displaying back-and-forth 
sideways swaying motion during their life time. For all 
the RBEs that have a duration of more than two time 
steps, we measured how many times there is a change in 
direction. We found that for both Ha and Ca II 8542 
RBEs, the change in direction happens on average 1.7 
times during their lifetime. 

Figure [5] shows the variation of the Doppler velocity 
and Doppler width for RBEs along their length as color- 
coded lines drawn at their spatial location in the FOV. 




Figure 5. Histograms for measurements of the transverse motion of RBEs: (a) maximum transverse displacement or distance between 
the two most extreme positions during the lifetime of an RBE, (c) maximum transverse velocity of all measured velocities during the RBE 
lifetime, and (d) all measured transverse velocities over all time steps. The scatter plot in panel (b) shows the total transverse distance 
covered by an RBE during its lifetime against the maximum transverse displacement of (a). The solid line marks the line where the two 
measures are equal, the dotted line marks where the total transverse displacement is 50% larger than the maximum transverse displacement. 



All detected RBEs are drawn in successive order, with 
RBEs detected early in the time series drawn first so 
that they are covered by later RBEs that occur at the 
same spatial location. The RBE occurrence in the target 
area is dominated by the magnetic region which gives a 
clear sense of the orientation of the RBEs and where they 
originate. This representation gives a clear visual impres- 
sion of many RBEs increasing in both Doppler velocity 
and Doppler width towards their top. This trend is most 
clear for the Ha RBEs. It should be noted that despite 
the visual trend, there is a significant number of RBEs 
that have a more erratic variation of Doppler properties 
along their axis and even examples of decreasing trends. 

There is a clear correlation between increased Doppler 
width for increasing Doppler velocity, this is illustrated 
in Fig. [3 

4.3. The RBE Ca II 8542-Ha Connection 

From the initial detection of 2717 and 2099 events in 
Ca II 8542 and Ha, respectively, a total of 1890 and 
1601 events combined to form 507 and 453 multi-frame 
RBEs. These chains from both spectral lines were then 
compared to one another in order to see if there is a 
connection between the Ca II 8542 and Ha RBEs. The 
same method that was used for linking separate events to 
multi-frame RBE chains was employed to link the chains 
in the two spectral lines. The automated method man- 



aged to find about 1/3 of the multi- frame RBEs in both 
spectral lines simultaneously resulting in 169 Ca II 8542- 
Ha RBE chains. Figure \Vp shows the distribution of 
these 169 chains, drawn on a wideband image. The 
Ha parts of the chains are drawn as green lines, the 
Ca II 8542 parts as red lines, and the region where they 
overlap is drawn in yellow. 

Our detection technique finds that a large fraction of 
the RBEs are detected and matched up in both lines, 
with a small fraction not matching up. However, from 
visual inspection of the dataset, it becomes clear that for 
virtually all RBEs detected in one spectral line a corre- 
sponding RBE signature is observed in the other spectral 
line. This is illustrated in Fig. and b, where close-ups 
of one time step are shown for both lines with the de- 
tected events in the other spectral line over-plotted. The 
automated detection method managed to link 7 RBEs 
in both spectral lines: 4 are found in both lines for this 
time step (marked as 1,2,3, and 4), 1 Ca II 8542 RBE 
is linked to an Ha RBE in another time step (marked 
'a'), and 2 Ha RBEs are linked to Ca II 8542 RBEs in 
another time step (marked 'b' and 'c'). The other de- 
tected RBEs were not flagged as connections (marked 
with small stars). Visual inspection however suggests a 
much higher connection rate which becomes even more 
apparent in the composite map in panel c. Here, Ha and 





Figure 6. Doppler properties of RBEs as a function of their position in the FOV, He* RBEs at left, Ca II 8542 RBEs at right. Color-coded 
Dopplcr velocities (top panels) and Doppler widths (bottom panels) drawn on background images of (a) Ha line centre, (b) Ca II 8542 line 
centre, (c) Ha wideband, and (d) Ca II 8542 —20 km s -1 detection map. 



Ca II 8542 detection maps are combined in a color im- 
age where bright Ha features are blue-green, and bright 
Ca II 8542 Doppler features are red. Regions where these 
features overlap (i.e., overlapping RBEs) become white. 
This map suggests that all RBEs have signal in both Ha 
and Ca II 8542. 

Panels (d), (e), and (f) of Fig.[8]offer a close-up centred 
on one of the longer RBEs. This long RBE was detected 
in both Ha and Ca II 8542 (see panels (a) and (b)), but 
the automated method did not connect them. The im- 
ages (and in particular the composite image) show that 
there is considerable overlap between the two lines. The 
other long but thinner RBE to the right is also detected 
in both lines but not flagged as connections - the com- 
posite image leaves no doubt that the two detections are 
of the same event. And the small RBE in the upper left 
was only detected for Ha for this time step (marked V, 
connected to a Ca II 8542 detection in the next time step) 
but there is indeed RBE signal present also in Ca II 8542. 



Figure[8]demonstrates that Ca II 8542 RBEs have their 
lower part closer to the magnetic regions, while Ha RBEs 
extend to larger distances away from these regions. This 
trend of Ca II 8542 RBEs being closer to the magnetic 
regions is also illustrated in Fig.[Tfc: most of the red lines 
(Ca II 8542) are at the base of the green Ha RBEs. The 
single RBE in Fig. [2] is also found in Ca II 8542 with the 
latter feature located closer to the magnetic bright points 
in the upper right (where the RBE started to appear a 
few time steps earlier) than in Ha. 

When comparing the statistics of RBEs connected in 
both spectral lines to the ones that are not connected, 
two of the properties stand out. The average lifetime 
of an RBE seen in both spectral lines is 83.9 s whereas 
in the RBEs only detected in one of the lines, the aver- 
age lifetime is 45.6 s and 46.1 s for Ca II 8542 and Ha, 
respectively. Also the average length of RBEs seen in 
both spectral lines, 3.1 Mm and 3.7 Mm for Ca II 8542 
and Ha, respectively, is larger than that of the RBEs 
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Figure 7. Scatter plot of average Doppler width vs. average 
Doppler velocity of all RBEs, represented as a 2-dimensional den- 
sity functions. Darker regions represent the densest areas with data 
points. 

only observed in one spectral line, 2.8 Mm and 3.4 Mm 
for Ca II 8542 and Ha, respectively. This suggests that a 
successful connection between the two lines is more likely 
for RBEs that are longer lived and have longer spatial 
extent. 

The statistics for Doppler velocities and Doppler 
widths for the connected RBEs shows the same trend 
as for RBEs detected in only one line: Ha RBEs have 
generally larger Doppler velocity and larger width. 

4.4. Scattered Events 

Figures Q] and |5] show that apart from the dense concen- 
trations of RBEs associated with regions of enhanced po- 
larization, there are a significant number of events scat- 
tered over the FOV. Here we present statistical properties 
of these isolated events by masking out the three regions 
with densest RBE concentrations. 

Figure O shows histograms for length, Doppler velocity 
and Doppler width. The histogram for RBE length peaks 
just below 3 Mm, just like for the multi-frame RBEs 
in Fig. 2K, but there is no extended tail towards longer 
lengths (note that some of the scattered RBEs show tem- 
poral evolution and are part of both Fig. 2] and EI- The 
histograms for Doppler velocity and width are similar as 
for multi-frame RBEs. 

Figure [TO] shows sample images of two scattered Ha 
RBEs. The top panels are from an event that appeared 
one time step earlier in the far blue wing in the middle 
of a granule. The RBE could be followed for 4 time 
steps so it has a life time of more than 36 s. In the 
left panel, the RBE appears to cross a granule, while in 
the images closer to line core (for example the middle 
panel), it seems to be part of a structure that is rooted 
in the intergranular lane to the right of the granule. No 
photospheric bright point or enhanced polarization signal 
that is obviously associated with this event is present in 
this area. 

The bottom panels show an RBE that was visible for 8 
time steps, and like the previous example seems to cross 
granules and is not obviously connected to magnetic re- 
gions in the vicinity. The line core image however sug- 
gests that the RBE is associated with the dense fibrils 
that are rooted in the extended magnetic region to the 
upper right that dominated the FOV in Fig [1] The clos- 



est photospheric bright points along the direction of the 
RBE main axis are more than 10" away. 

5. DISCUSSION 

We studied the disk-counterparts of type II spicules 
that are observed as short-lived, elongated features in 
the blue wings of strong absorption lines like the Ha and 
Ca II 8542 spectral lines. These features are also known 
as rapid blue-shifted excursions due to their asymmet- 
ric, blue-ward spectral profile. The work presented here 
is in many ways a continuation and elaboration of the 
work of Paper I. We observed a similar target, a coro- 
nal hole, and employed a program of simultaneous Ha 
and Ca II 8542 spectral scans. The observation program 
was specifically tailored for the observation of RBEs with 
reasonably high temporal cadence and sufficient resolu- 
tion for robust spectral analysis. This allowed us to more 
directly compare the two spectral diagnostics, a compar- 
ison which could only be done indirectly from the data 
employed in Paper I, where RBEs were studied in two 
separate but sequentially recorded single-line time series 
in Ha and Ca II 8542. Furthermore, we have extended 
the method of automated detection of RBEs to allow 
for a temporal connection between RBE detections so 
that properties such as lifetimes and transversal motion 
could be studied more rigorously on an extended statis- 
tical sample. In Paper I, these properties were measured 
manually on a limited subsample of the RBE detections. 

5.1. Statistical properties of RBEs. 

We present histograms on RBE lengths, Doppler veloc- 
ities and Doppler widths for 1890 Ca II 8542 and 1601 
Ha RBEs in Fig. [4] These measurements agree very well 
with those presented in Paper I. We confirm that a larger 
number of longer RBEs (>3.5 Mm) can be observed in 
Ha, while the shortest RBEs (<2 Mm) are observed in 
Ca II 8542. We confirm that in Ha, one observes larger 
Doppler velocities and larger Doppler widths than for 
Ca II 8542 RBEs. 

These RBEs were connected in temporal chains and 
combined to form 507 Ca II 8542 and 453 Ha multi- 
frame RBEs for which we could determine the life time 
and transversal motion. This is an increase in statistical 
sample size of more than an order of magnitude for this 
type of measurements. 

Pereira et al. (private communication) measure statis- 
tical properties of 111 type II spicules from two Hinode 
Ca II H limb time series. Their results allow for a more 
detailed comparison between RBE properties and spicule 
type II properties than was possible for Paper I. They 
find life times ranging between 20-160 s, with the peak 
of the distribution around 1 min and an average of about 
80 s. For our sample of RBEs, we find a similar range 
in life times but the peak of the distribution well below 
1 min and an average of 50 s. We note that for the RBEs 
that can be identified in both Ca II 8542 and Ha, we find 
an average life time that agrees well with the Pereira et 
al. result. 

Pereira et al. find for the transverse motion the av- 
erage velocity ranging between and 20 km s" 1 , with 
the peak of the distribution at about 6 km s _1 . We find 
a similar range in velocities, but most of the measure- 
ments for RBEs yield lower velocities, with the average 
at -4 km s _1 for Ca II 8542 and -5 km s _1 for Ha. The 
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Figure 8. Examples of RBE detections, (a): Ha —45 km s 1 map used for RBE detection, RBEs are white. Red lines mark Ca II 8542 
detections from the same time step, (b): Ca II 8542 —20 km s —1 detection map, RBEs are white. Green lines mark Ha detections from the 
same time step, (c): composite image of the Ha and Ca II 8542 detection maps: the Ca II 8542 image filled the red channel, the Ha image 
filled the blue and green channels of the true color RGB image. Detections marked 1-4 are RBEs connected in both Ha and Ca II 8542 for 
this time step, detections marked a-c are RBEs connected to the other spectral line in a different time step, detections marked with stars 
are RBEs that are not identified in the other line. The three bottom right panels show a closer view of one RBE marked with a larger star 
at (x,y) = (39, 28): (d) the Ha detection map, (e) the Ca II 8542 detection map, and (f) the Ha-Ca II 8542 composite image. 

measurements of the transverse motion of 35 RBEs in Pa- 
per I, resulted in a transverse velocity of order 8 km s _1 . 
These Ha RBEs were selected manually but we con- 
firmed that result by applying the method of automated 
measurement of transverse motion to the 2008 data set 
and found an average of 7.6 km s _1 for 179 Ha RBEs. 
We attribute this higher value for the transverse velocity 
to the higher cadence of the 2008 Ha data set which was 
almost two times faster than the 2010 data. This allows 
for the measurement of higher velocity motions because 
the upper limit of maximum allowed distance between 
time steps is larger. For the 2008 Ca II 8542 data set, 
which is at similar cadence as the 2010 data, we find an 
average of 4.4 km s" 1 for 94 RBEs. For the 03-July-2010 
data set, we find average transverse velocities that are 
very similar as for 27-Jun-2010. 



iDe Pontieu et al.l (|2007| ) discovered Alfvenic waves in 
the chromosphere using type II spicules as tracers for os- 
cillatory motions. The velocity amplitude was measured 
for 94 spicules which can be compared to the distribu- 
tion of the peak transverse velocity in Fig. [SJ:. They 
found a range between and 30 km s _1 , with the peak 
of the distribution around 15 km s _1 . Like for the aver- 
age transverse velocity, RBEs display a comparable range 
but mostly lower values for the peak transverse velocity 
as compared to type II spicules at the limb. Continuing 
on the arguments of the previous paragraph, we note that 
the measurements of the RBE transverse velocities seem 
to be limited by the cadence of the time series. Faster 
cadence observations are needed to properly resolve the 
distribution of transverse velocities of RBEs. 

For the transverse displacement, Pereira et al. find 
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Figure 9. Histograms for (a) length, (b) Doppler velocity, and 
(c) Doppler width for RBEs scattered over the FOV, i.e., not as- 
sociated with enhanced polarization signal (see Fig[T]i). 

values ranging between and 1.4 Mm, with most values 
between 0.1 and 0.5 M m, in close agreement with the 
iDe Pontieu et al.l (|2007l ) measurements. We find a total 
transverse displacement of less then 1 Mm for RBEs, 
with most RBEs covering a distance of less then 0.2 Mm 
during their life time. 

Even though the shapes of the statistical distributions 
of properties of RBEs and type II spicules observed at 
the limb differ in detail, the similarities between the dis- 
tributions provide sufficient basis for the interpretation 
of RBEs being the disk counterparts of type II spicules, 
as already put forward in Paper I. In addition to similar- 
ity in statistics, we observe that RBEs rapidly disappear 
or fade, and that RBEs show both swaying motion and 
apparent motion away from magnetic regions that can 
be interpreted as an upward motion. There are no signs 
of flows returning to the surface associated with RBEs. 
This kind of dynamical behavior is one of the defining 
characteristics of type II spicules observed in Ca II H at 
the limb. 

The most striking difference between RBEs and type 
II spicules is the distribution of lengths. Pereira et al. 



measure a range between 1.5 and 11.5 Mm, with the av- 
erage length around 6 Mm. We find for RBEs an average 
length around 3 Mm and only few Ha RBEs measure 
up to 6 Mm. This difference can be attributed to the 
fact that above the limb, faint signals towards the top 
of spicules can still be reliably detected against the dark 
background. Against the disk however, these top parts of 
spicules are likely to have too little opacity to be detected 
as an integral part of RBEs. So there is likely an intrinsic 
bias for shorter lengths measured for RBEs. For length 
measurements of type II spicules at the limb, there is ar- 
guably also a selection effect for detecting longer spicules 
that stick out above the dense forest of spicules encoun- 
tered in limb observations. The same projection effect 
could possibly explain the trend of higher transverse ve- 
locity and transverse displacement measured for spicules 
at the limb: tall spicules that extend to greater heights 
move in a lower density environment, so for equal wave 
energy flux they can be expected to display higher veloc- 
ities and thus cover greater transverse distances. 

5.2. RBE detection rate. 

In Paper I, in the 40 min Ca II 8542 time sequence with 
cadence 11 s, 413 RBEs were detected from Dopplermaps 
at 30 km s _1 . With the 215 time steps, this corresponds 
to a detection rate of 1.9 RBEs per time step. In the 
24 min Ha time sequence with cadence 6.7 s, 608 RBEs 
were detected from difference maps at 60 km s _1 . With 
the 245 time steps, this corresponds to a detection rate of 
2.5 Ha RBEs per time step. We note that detection ve- 
locities of 30 km s" 1 for Ca II 8542 and 60 km s" 1 for Ha 
are on the extreme end of the average Doppler velocities 
found for RBEs, see Fig. [4] (or Fig. 10 in Paper I). Given 
the typical values for the Doppler width and velocity it is 
likely that at these extreme detection velocities, a signif- 
icant fraction of RBEs have absorption signals that are 
too faint to be detected by the automated method. In 
this work we investigated the detection rate as a function 
of the detection velocity used in the automated method. 
We find similar detection rates as Paper I at the detec- 
tion velocities used in Paper I: 3.4 RBEs per time step for 
Ca II 8542 and 1.8 RBEs per time step for Ha. Running 
the automated detection method on maps constructed at 
lower Doppler velocity, we find a significant increase in 
the detection rate, see Fig. [3] At 30 km s" 1 we had a 
detection rate of 48 Ha RBEs per time step, an increase 
in the number of detected events of more than a factor 
25. For Ca II 8542 we see a similar increase in the de- 
tection rate towards lower Doppler velocity but we found 
that at the lowest Doppler velocity possible in our data 
set, 12 km s -1 , the fraction of false detections was too 
large and the detections are not reliable. For the detec- 
tion velocity for which we present RBE statistics, 20 and 
45 km s _1 for Ca II 8542 and Ha, we see an increase in 
the detection rate with a factor of about 6, to 14.2 and 
11.1 RBEs per time step, as compared to the Paper I de- 
tection velocity. This resulted in a significant increase of 
the RBE sample used for detailed analysis in this work. 

We find similar trends of increased RBE detections at 
lower Doppler velocity for the 03-Jul-2010 data set and 
the 15-Jun-2008 data sets used in Paper I. 

We also note that the automated detection method 
does not identify all RBEs present in the data. We es- 
timate that for the 45 km s _1 Ha detection maps, only 
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Figure 10. Sample Ha images of two scattered RBEs. Left panels are images in the blue wing at 
and right panels in line centre. 
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about half of the events that are visually identified as 
genuine RBEs are registered by the automated detec- 
tion. Running the detection routine on lower Doppler 
velocity maps results in a lower relative number of unde- 
tected RBEs (about a quarter remains undetected) but 
we must underline that we considerably underestimate 
t he abundance of R BEs in our data. 

Judg e et "ail (|2011h proposed that type II spicules cor- 
respond to warps in two-dimensional sheet structures 
as an alternative to the main-stream interpretation of 
spicules being tube-like structures. One of their ar- 
guments in favor of sheet-like structures is that RBEs 
(which are hard to interpret as " sheets" ) are not the disk 
counterparts of type II spicules because the occurrence 
rate of RBEs is at least an order of magnitude lower than 
type II spicules observed at the limb. From the number 
of detected RBEs in Paper I, they extrapolate to an esti- 
mate of 10 5 RBEs on the Sun at any given time. This is 
then orders of mag nitudes lower than estim ates for type 
II spic ules, 2 x 10 7 pudge fc Carlssonll2010D . and lBeckersI 
(1968) 's estimate of 10 6 for what is now considered as 
classical spicules. 

The detection rate of RBEs at lower Doppler velocity 
we find here agrees with the estimate of a total number 
of ~ 10 5 RBEs (see Fig. [3j). However, as argued above, 
the RBE detection rate we find here is an absolute lower 
limit that significantly underestimates the actual number 
of spicules for several reasons. 

First, the automated detection rate rejects a significant 
fraction of real RBEs to avoid false positives. 

Secondly, there are clear indications that many spicules 
lack opacity to show a disk counterpart signal. Detecting 
an RBE-like feature against the disk requires observa- 
tions of excellent quality, and most importantly requires 
that the RBE has a certain degree of enhanced opacity. 
Our disk observations indicate that the opacity of RBEs 
varies strongly between the Ca II 8542 and Ha lines, 



with the former showing significantly fewer and shorter 
features (concentrated towards the bottom of the Ha fea- 
tures). This is compatible with an ever decreasing opac- 
ity in chromospheric lines as one moves away from the 
footpoints. This is not surprising because limb spicules 
also show a very strong decrease in intensity (with scale 
height of order 2000 km) over their full length. It is 
thus highly likely that there is a significant number of 
faint RBEs that cannot be detected even in Ha. We 
note that RBEs are very difficult to observe: detecting 
RBEs was impossible before the advent of subarcsecond, 
high cadence, high signal-to-noise imaging spectroscopy 
with Fabry-Perot interferometers. It is quite possible 
that even higher quality observations would increase the 
detection rate significantly. 

Thirdly, RBE detection is subject to strong velocity 
filtering: due to the unknown angle between the velocity 
vector and line of sight vector, our method always misses 
features that are significantly inclined from the line of 
sight. This is especially the case because we cannot dis- 
tinguish RBEs from other chromospheric features at ve- 
locities closer to the core. The increase of the number of 
RBE detections for lower Doppler velocity as illustrated 
in Fig. 3 indicates that there might be a significant num- 
ber of RBEs at even lower Doppler velocity that are not 
detected by our method. Detection of type II spicules at 
the limb suffer much less from these latter two effects: 
it is relatively easy to detect faint signals against the 
off-limb dark background (with space-based instruments 
like Hinode that do not suffer from seeing effects intro- 
duced by the Earth's atmosphere) and velocity filtering 
is much less an issue for the relatively wide pass-band of 
the Hinode/SOT Ca H filter. 

In the discussion of Paper I, it was argued that the 
detection rate of RBEs on the disk correspond to a de- 
tection of 1.9 RBEs per linear arcsec if they were ob- 
served at the limb. This was considered as being in 
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reasonable agreement with an order of magnitude esti- 
mate of the number of type II spicu les per linear arc- 
sec from the iDe Pontieu et all (|2007f ) observations. In 
this study, we confirm the results from Paper I and find 
that the detection rate extrapolates to ~ 1 — 2 x 10 5 
RBEs over the whole Sun. We caution that all of these 
occurrence rate estimates, in c luding those of iBeckersj 
(1968) and Uudee fc Carlssonl (I201OT are based on in- 
complete datasets, limiting assumptions, observational 
limitations, and should thus be take n with a grain of 
salt. F or example, it is not clear from lJudge fe Carlssonl 
(2010) how deterministic their occurrence rate of 10 7 
type II spicules really is. The visual comparison in their 
paper is unconvincing with their simulations showing 
many more spicules than observations do: ~ 10 7 type 
II spicules would correspond to more than 100 spicules 
per linear arcsec at the limb or more than 5 spicules 
per Hinode/SOT pixel, in contradiction to observations. 
This implies that t he 10 7 type II spicules estimated by 
Uudee fc Carlssor] (pOlOf ). is likely a significant overesti- 
mation of the number of type II spicules at the limb. In 
any case, their argumentation for that number is based 
on a series of assumptions and models whose effect on the 
estimated number of spicules has not been studied or es- 
timated. This sheds doubt on the robustness of their 
determined detection rate. 

We conclude that the rejection of lJudge et al.l (|2011D of 
RBEs being the disk counterparts of type II spicules on 
the basis of too low occurrence rate is unjustified. There 
may well be sheet-like features in the atmosphere, but 
challenging the "tube hypothesis" based on rejection of 
the RBEs as disk counterparts of spicules is not convinc- 
ing for the reasons outlined above. 

5.3. Acceleration along RBEs. 

Paper I reported a systematic variation of the Doppler 
velocities and widths: many RBEs were found to show 
an increase of Doppler shifts and widths from the foot- 
point to the top. We confirm this trend in the 2010 data 
which is visually evident from Fig. [5] which shows the 
Doppler measurements as a function of their position in 
the FOV. In this representation, the orientation of (most 
of) the RBEs can be comprehended with respect to their 
footpoints which appear to be rooted in the magnetic 
field concentrations. The trend of increasing Doppler ve- 
locity and width for many RBEs is particularly striking 
for the Ha RBEs. We note that there is also a significant 
number of RBEs displaying a more erratic variation of 
their Doppler properties. 

Irrespective of the variation along the RBE main axis, 
we report a more general trend of increasing Doppler 
width for Doppler velocity in Fig. [7] This trend illus- 
trates the characteristic spectral profile at the location 
of RBEs which consists of a wide blue-ward asymme- 
try (from which the RBE profile is then extracted by 
subtraction of a background profile). This is different 
from profiles with a narrow, isolated absorption compo- 
nent which can for ex ample sometimes be observed for 
on-disk coronal rain (lAntolin fc Rouppe van der Voortl 
|20H lAntolin et all 120 111 ). Sometimes coronal rain con- 
densations appear as separated absorption profiles in the 
far wings of Ha. This is unusual for RBEs, even for 
RBEs that are oriented along the line-of-sight and ap- 
pear as small, roundish features with strong absorption 



in the blue wing (dubbed "black beads" in Paper I), the 
spectral profile displays enhanced absorption throughout 
the blue extending to large Doppler offset. 

In any case, the correlation between Doppler velocity 
and widths and tendency towards increased values for 
both parameters towards the top of the RBEs may sug- 
gest a link between the acceleration and heating mecha- 
nism involved in RBEs, if we speculate that the Doppler 
widths are related to temperature. Whatever the cause of 
this intriguing correlation, it provides strict constraints 
for any theoretical model for spicules. 

5.4. Comparing Ha and Ca II 854-2 RBEs. 

From the sequentially recorded Ha and Ca II 8542 
data sets of Paper I, it was observed that Ca II 8542 
RBEs are located closer to magnetic regions and it was 
inferred that the Ca II 8542 line samples the lower part 
of spicules. In our data sets with co-temporal Ha and 
Ca II 8542 observations, we can directly confirm that 
the two lines sample different parts of RBEs, with the 
Ca II 8542 line sampling the bottom part of the spicule, 
and Ha the upper part, and both lines sampling an over- 
lapping region in the middle (see Fig|8]). Even though the 
automated detection method does not nearly identify all 
RBEs in both lines, we find from visual inspection that 
virtually all RBEs display signal in both lines. 

This observation strengthens the explanation for the 
differences in Doppler velocity measured in Ha and 
Ca II 8542 (see Fig. [4]) as put forward in Paper I: Ha 
is sampling the top part of spicules where lower density 
plasma is propelled to higher velocities. 

The observation of the Ca II 8542 and Ha lines sam- 
pling different parts of RBEs also removes some of the 
discrepancy between length measurements of RBEs and 
type II spicules at the limb: we can infer that the ac- 
tual length of the RBEs is longer than measured in the 
two lines separately. In addition, of course, it remains 
the case that we cannot expect to completely remove the 
length discrepancy between RBEs and type II spicules. 
This is because off-limb measurements allow for more re- 
liable measurements of faint signals. 

5.5. Scattered RBEs. 

We find most RBEs to be concentrated near regions 
with enhanced magnetic fields. In addition to these 
dense concentrations, we find a number of more isolated 
RBEs that are scattered around the FOV. When com- 
paring statistical properties of length, Doppler velocity 
and Doppler width, we find that these isolated events are 
not particularly different from the RBEs in denser re- 
gions, except that we do not find a population of longer 
RBEs. For these isolated RBEs, it is often unclear if they 
have any association with regions of enhanced magnetic 
fields, although we note that our diagnostic, Ca II 8542 
—600 mA Stokes V maps, is not a particularly sensitive 
measure for magnetic fields. One might speculate that 
a more sensitive diagnostic might reveal a connection to 
weaker and possibly more extended, magnetic fields. 

From fixed wavelength, —1300 mA Ha observations 
from the 1.6 m Ne w Solar Telescope ( NST) in Big Bear 
Solar Observatory, iGoode et all (|20101 ) report the obser- 
vation of small jet-like features that appear to originate 
from intergranular lanes. The authors claim that these 
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features are different from RBEs, as they are not un- 
equivocally tied to strong magnetic field concentrations. 

Figure [TOl shows blue-wing Ha images of two examples 
of scattered RBEs. The left panels show far blue wing im- 
ages of these events. In these images, these events seem 
to match the descriptio n of the jet-like features described 
bv lGoode et~aH (I20T0I) : narrow features that cross gran- 
ules and are not obviously connected to magnetic regions 
in the vicinity. The spectro-temporal information that 
we have for this data however reveal that these events 
are RBEs. Even though the bottom example of Fig. [TO] 
is not directly associated with photospheric bright points, 
it seems plausible that this feature is related to the ex- 
tended magnetic region in the upper-right part of the 
F OV. 

iGoode et al.l (|2010l ) report a typical length of 1 Mm 
and a typical width of 0.2 Mm for the jet-like features in 
their blue-wing Ha data. The width agrees well what we 
find for RBEs, the length is compatible but on the short 
side. It should be noted that the length in IGoode et al.l 
(2010) is measured at a fixed wavelength position at 
60 km s _1 Doppler offset, while we measure RBE lengths 
taking the spectral profile into account. This naturally 
leads to a trend of measuring longer lengths. We con- 
firmed visually that RBE lengths at fixed wavelength 
filtergrams have a trend of becoming shorter at higher 
Doppler offset. 

Given the striking sim ilarities, we i nterpr et the jet- 
like features reported by IGoode et a l. (2010) as similar 
features to what we observe as isolated RBEs scattered 
over the FOV. 



6. CONCLUDING REMARKS 

The statistical properties resulting from this dataset, 
along with better estimated occurrence rates, further 
strengthen the interpretation of RBEs as the disk- 
counterpart of type II spicules as first proposed by 
Langang en et al.1 (|2008) and further established in Pa- 
per I. Given their important role in understanding 
the outer solar atmosphere, either as seemingly pas- 
sive tracer for the perva sive presence of Alfvenic waves 
(|De Pontieu et al.l l2007h or more directly as a poten- 
tial source for mass loa ding and heating of the corona 
(|De Pontieu et al.l 1201 If ), further observational study of 
type II spicules is highly desired in order to advance 
knowledge of their physical nature. An essential property 
that is not part of this work but can readily be addressed 
from the analysis method employed on this time series 
is the spectral evolution of RBEs during their lifetime. 
Such detailed study of the temporal evolution of RBEs 
is currently under way and will be the subject of a forth- 
coming paper. We further note that observations from 
the Solar Dynamics Observatory of this target provide 
the possibility to investigate the response of transition 



region and coronal diagnostics to RBEs. 
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